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Q_i' A universal seesaw mass matrix model, which gives successful description of 

Oh: 



quark mass matrix in terms of lepton masses, yields three "sterile" neutrinos 
vf, which compose pseudo-Dirac neutrinos ~ {vi it f|)/\/2 together with 



^ ■ the active neutrinos Vi (i = e,fi,T). The solar and atmospheric neutrino 



data are explained by the mixings i/^ <-> f| and <-> f^, respectively. In 

spite of such observations of the large mixing sin'^ 20 ~ 1 in the disappearance 

experiments, effective mixing parameters sin'^ 26°^^ in appearance experiments 

i^a 1^(3 (o, /3 = e, /X, r) are highly suppressed. 
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It has, for many years, been a longed-for goal in the particle physics to give a unified 
description of masses and mixings of quarks and leptons. The most crucial clue to the unified 
description is in the investigation of the neutrino mass matrix. Recent remarkable progress 
in the experimental studies of the solar, atmospheric and terrestrial neutrinos has put the 
realistic investigation of the neutrino mass matrix model within our reach. 

As a promising model which gives a unified description of quark and lepton mass matrices, 
the so-called "universal seesaw" mass matrix model has recently revived. The model has 
hypothetical fermions Fi = Ui, Di, Ni and Ei {i = 1,2,3) in addition to the conventional 
quarks and leptons fi = Ui, di, Ui, Cj, where these fermions belong to fi = (2, 1), //? = (1, 2), 
Fl = (1, 1) and Fr = (1, 1) of SU(2)ixSU(2)/j. Note that the neutral leptons Ni are the 
so-called "sterile" fermions because they are singlets of SU (2)^x811(2)/? and do not have 
U(l)-charges. It has recently been pointed out by one of the authors (YK) that such a 
model can yield three massless sterile neutrinos u- = {Nu — Nr^)/\/2 [Nr = {NrY = CiV^] 
in a limit, and the sterile neutrinos compose pseudo-Dirac neutrinos uf^ ^ (z/j±z/j'*)/-\/2 
with the masses m(i/f_^) ~ m(z/f_!) together with the conventional left-handed neutrinos 

The purposes of the present paper are to investigate an explicit model with three sterile 
neutrinos within the framework of the universal seesaw mass matrix model and to give 
numerical predictions for neutrino masses and mixings. The solar 0] and atmospheric 
data will be explained by the mixings Ve ^ vl and ^ w^, respectively. A similar 
idea has recently been proposed by Geiser [^], who has introduced three sterile neutrinos 
phenomenologically. However, in the present model, our sterile neutrinos u- come from 
the sterile fermions Ni in the universal seesaw mass matrix model, so that their masses and 
mixings are strongly constrained by the parameters determined at the charged lepton masses 
and quark masses. 

The seesaw mechanism was first proposed [j^ in order to answer the question of why 
neutrino masses are so invisibly small. Then, in order to understand that the observed quark 
and lepton masses are considerably smaller than the electroweak scale A^, the mechanism 
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was applied to the quarks [jT| : A would-be seesaw mass matrix for (/, F) is expressed as 



M 



rriL 




( 






= mo 















(1) 



where the matrices Zl, Zr and Yf are of the order one. The matrices itil and nipi take 
universal structures for quarks and leptons. Only the heavy fermion matrix Mp takes a 
structure dependent on f = u, d, z/, e. For the case A ^ k ^ 1, the mass matrix (1) leads to 
the well-known seesaw expression of the 3x3 mass matrix for the fermions /: 



Mf ~ -mLMp^rriR . 



(2) 



However, the observation of the top quark of 1994 ^ once aroused a doubt that the observed 
fact rrit ~ 180 GeV ~ = 0{mL) was not able to be understood from the universal seesaw 
mass matrix scenario, because nit ~ 0{'mi) means Mp^rriR ~ 0{1). For this question, a 
recent study @,|l3] has given the answer that the universal seesaw scenario is still applicable 
to top quark if we put an additional constraint 



detMp = , 



(3) 



on the up-quark sector {F = U). The light quarks {qLiARi) acquire the well-known seesaw 
masses of the order of A^Aij/Ag though the heavy fermion masses m{Qpi, Qui) ~ A^- [A^ = 
mo = O^itil), Afi = kuiq = 0{mji), and As = Xtjiq = 0{Mp)], while the third up-quark 
ul3 acquires a mass of the order of Aj;^ together with the partner Urs since the seesaw 
mechanism does not work for the fermions (ulsjUrs) because of m^UpsjUfis) = 0. Thus, 
by regarding the fermion state (mls, Urs) [not ('Ul3,m_r3)] as the top quark, we will easily be 
able to understand why only top quark t acquires the mass m^ ~ 0{mL). A suitable choice 
of the matrix forms of Zp, Zr and Yf can give reasonable values of the quark masses 
and Cabibbo-Kobayashi-Maskawa [|1T| (CKM) matrix parameters in terms of charged lepton 
masses. 

On the other hand, in the conventional universal seesaw mass matrix model [Q, the light 
neutrino mass matrix Mjj is given by [|12 
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(4) 



so that the smallness of the neutrino masses rrii, is understood by a large value of k = 
Ar/Al {k > 1) because of m,y ~ A\/ As ~ {l / K)mchaTged lepton- The case (4) can be 
brought, for example, by assuming that the neutral leptons and N^i which are singlets 
of SU(2)2,xSU(2)ii; and have no U(l) charges acquire Majorana masses Mm ~ in addition 



to the Dirac masses Mo = Mn ~ As and assuming Mm = Mc ||T3 

( . . . \ 



mi 
m]J 
tur Mm Md 



m 



T M, 



(5) 



where Mj, is the mass matrix which is sandwiched between {v i, 77^, N^, N'j^) and 
{u'1,ur, Nl, Nr)'^ , and z/£ = {i'lY = Cv\ and so on. The model given in Ref. [0 can 
readily lead to a large mixing between and Vr [or z/g and z/^] with m\ <^ m!^ — rrig [or 
m\ ~ 1712 ^ "^3] ) so that the model is favorable to the explanation of the atmospheric neu- 
trino data [H which have suggested a large neutrino mixing ^ v^- However, the model is 
hard to give a simultaneous explanation of the atmospheric ||^ and solar |^ neutrino data. 

We want to give a simultaneous explanation of the atmospheric and solar neutrino data. 
The essential idea is as follows [0: In the present model, there is no quantum number which 
distinguishes A''^ from Nr^ because both fields are SU(2)ixSU(2)/j singlets and do not have 
U(l) charges. Therefore, if vj^ (ur) acquire masses rriL {niR) together with the partners Nr 
(Nl), they may also acquire masses m'^^ (m'^) together with the partners [Nf^). Then, 
the mass matrix for the neutrino sector is given by 



M 



^ m'l rriL ^ 
m'^ 
m'l niR Mm Mr, 



ml m'j^ Ml Mm 



(6) 



By rotating the fields (z/£, vr, N^, Nr) by 



R. 



■34 



^10 ^ 
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the 12 X 12 mass matrix (6) is expressed as 

/ 



M' = RziMRl^ 



V2AJ 




V 
V2AI V2Ar Mm - Md 
V2ml V2mn 



rriL ^ 



m 



■T 



Mm + Md 



where we have used M^ = Md, and 



2Al = m'r - iTLL , 2Ar = mR- m'l 



R ' 



2mL = tn'i^ + rriL , 2mR = rriR + . 

In the Nl ^ Nr symmetric hmit, i.e., in the hmit of m'^^ = rriL, itl'r = rriR and 
Md = Mjv, the mass matrices for the neutrinos u^, and ur are given by 

M{u') = , 



M(N') ~ 2M, 



N 



where z/* and A^'' are defined by 

Furthermore, for a model with the relation tur oc rriL such as a model given in Ref. 
obtain 



MM = . 



(15) 



Only when we assume a sizable difference between m'^ and uil (and also between m'^ and 
ttir), we can obtain visible neutrino masses m{uL) 7^ and m^u^) ^ 0. 
For Ai, 7^ and A/j 7^ 0, by using a unitary matrix U12 
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[CC^ + SS^ = 1], we can express the mass matrix (8) as 

/ 



M" = Ui2M'U^^ 



A 

5 G 

B'^ Q 

G'^ 2Mjv 



\ 



where we have put Mm — Mr, = and 



A = V2{CAl - SAl) , 



(16) 



(17) 



(18) 



G = V2{S^mL + C^ml) , 



(19) 
(20) 



CS-^ = m^i(mL) 



-1 



(21) 



Since M^ ~ A5, G ~ Ar, 0{B) <^ A^j and 0{A) <^ A^, by using the seesaw approximation, 
we obtain the 9x9 mass matrix for states {ul, ur, u^) 



( 



A 

-G(2M^)-iG^ B 
A^ B^ 



(22) 
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together with the 3x3 mass matrix (12) for the neutral leptons A^*. Furthermore, for 
0{GM^^G^) = Aji/As = iK^/X)mo > 0{A),0{B), we obtain the 6 x 6 mass matrix for 

f A \ 



together with the 3x3 mass matrix for the right-handed neutrinos upt 



(23) 



M{ur) ~ -G(2Mjv)~'G' . (24) 

If we choose 0{A) ^ 0{B{G'^)~^MnG~^B'^), we obtain three sets of the pseudo-Dirac 
neutrinos (z/f_^, ufl) {i = e, /i,r) which are large mixing states between z/j and u-, i.e., 

^s^^(^.±^r), (25) 

and whose masses are almost degenerate. 

So far, we have not assumed explicit structures of the matrices Zl, and Yf. Here, in 
order to give a realistic numerical example, we adopt a model with special forms of Z^, Zji 
and Yf which can lead to successful quark masses and mixings. The model is based on the 
following working hypotheses |^: 

(i) The matrices Zl and Zr, which are universal for quarks and leptons, have the same 
structure: 

Zl = Zr = Z = diag(zi, Zs, z^) , (26) 

with z\ + zl = 1, where, for convenience, we have taken a basis on which the matrix Z 
is diagonal. 

(ii) The matrices Y/, which have structures dependent on the fermion sectors f = u,d, u, e, 
take a simple form [(unit matrix) + (a rank one matrix)]: 

Yf = l + 3bfX . (27) 

(iii) The rank one matrix X is given by a democratic form 
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1 1 1 



(28) 



on the family-basis where the matrix Z is diagonaL 

(iv) In order to fix the parameters 2;,, we tentatively take fog = for the charged lepton 
sector, so that the parameters Zi are given by 

1 (29) 



^Jrnl + + nir ' 

from Me ~ rriLMj^^mji = (K/\)mQZ ■ 1 ■ Z. 

The mass spectra are essentially characterized by the parameter bf. We take bu = 
— 1/3 for up-quark sector, because, at b^ = —1/3, we can obtain the maximal top-quark 
mass enhancement ^ (l/-\/3)wio, and a successful relation mu/mc — (3/4)(me/m^) 
independently of the value of k/X. The value of k/X is determined from the observed ratio 
rric/mt as k/A = 0.0198. Considering the successful relation mdims/rn^ ~ Ameim^/ml for 
bd — —1, we seek for the best fit point of bd = — e*^"^ <^ 1). The observed ratio 
iTLd/ms fixes the value (3d as (3d = 18°. Then we can obtain the reasonable quark mass 
ratios, not only m^/m^, mf/m'j, but also m^/mj: rriu = 0.00234 GeV, rric = 0.610 GeV, 
mt = 0.181 GeV, rud = 0.00475 GeV, = 0.0923 GeV, rub = 3.01 GeV, where we have 
taken (moK/A)^/ (mo/t/A)e = 3.02 in order to fit the observed quark mass values at /i = mz 
| 1^ . We also obtain the reasonable values of the CKM matrix parameters: |V^s| = 0.220, 



\Vcb\ = 0.0598, iKfel = 0.00330, \Vtd\ = 0.0155. (The value of |Kfe| is somewhat larger than 
the observed value. For the improvement of the numerical value, see Ref. |p.5|.) 

Stimulated by this phenomenological success in the present model for the quark masses 
and mixings, we apply the model with (26) - (28) to the neutrino mass matrix (6). For 
simplicity, we consider the case 



m 



X = (1 - £L)mL , m'^ = (1 - £R)mR . (30) 



Then, we obtain 



(31) 



(32) 



(33) 



Mn = moA(l + 3buX) . 



(34) 



What is of great interest to us is to evaluate the masses of the pseudo-Dirac neutrino states, 
whose mass matrix is given by (23): 



M(6) ^ —rnoisR + SL) 



^ z\ 



Z pY, 



(35) 



where 



p ~ 



1 e 



R 



-A . 



(36) 



V2eR + SL 

For Z » \pYu\, i.e., » |p(l + K)\, we obtain the masses m{uf^) of the pseudo-Dirac 



neutrinos uf^ 



'i±) - -^moisR + El) (zi ± P^^^ 



(37) 



Here, the double sign, ±, in the right-hand side of (37) does not always correspond to 
the double sign of uf^ in the left-hand side. For the mass eigenstates uf^ defined by (25), 
the order of the magnitudes of the mass eigenvalues are m{i^l) < m{i^^) < m(i/^+) < 
m(z/^!_) < m(z/^l) < m(z/^^) (see Fig. 1). On the other hand, the squared mass differences 
Am^(tf*) = \m'^{uf^) —m'^{i'fl)\ are given by 



Am {I'D ':^mQ{eR + El) \p\zi 



l + K)+ph^ — 



{zf - z]){z1 - zl 

where (i, j, k) is a cychc permutation of {i,j,k) = (1, 2, 3). 
As the numerical input, we use 



, (38) 



moj^^ = (mo^) = TrMe = (m, + m^ + m,)^=m^ = 1.8499 GeV , (39) 

(/c/A), = {k/\\ = 0.02 , (40) 

so that 

mo = 0.925 X lO" eV . (41) 

For simplicity, we seek for the solutions for the case = = e, so that the parameters 
in the neutrino masses and mixings are three, i.e., e, hp and A. We tentatively choose 
m(z/^4) ^ 5 eV, which is suggested from a cosmological model with cold+hot dark matter 
(CHDM) |jl6|. Then, for hy ~ —1, the parameter value e is given from (37) as 

£ = £^ = ~ 3.8 X 10"^^ . (42) 



The recent atmospheric data from Super-Kamiokande have suggested the neutrino os- 
cillation — > with Am^ ~ 2.2 x 10^'^ eV^ and sin^ 26* ~ 1. Since the value of Am^ 
in the case of z/^ — > z/^ is given by Am^(z/^ — z/^) ^ 2 x Am^(z/^ z/^) [|I8|, we put 
/S.m'^iy^^^) = Aml^^ ~ 5 x 10"^ eV^. Then, from (38), for b^, ~ —1, we obtain 

p - 0.057 , A ~ 4.3 X 10^ {k ~ 8.5 x 10^) . (43) 

The numerical results of the masses m(z/f*) = (m(z/f^) + m(z/ff ))/2 for the case of bi, ~ —1 
are as follows: 

m(z/f ) ~ 0.078 eV , m{uP') ^ 1.2 eV , m(z/f ) ~ 4.8 eV . (44) 

The value of Am^(z/^*) is highly sensitive to the parameter by. The large-angle Mikheyev- 
Smirnov-Wolfenstein (MSW) solutions [|1^] of solar neutrino problem, (Am^, sin^ 2^^) ~ 
(10~^ eV^,0.6), is ruled out for the case Ue — > ^'l 0]- Therefore, we use the vacuum solu- 
tion (Am^, sin^ 29) ~ (0.6 - 0.8 x 10"^° eV^ 1) [0,0 in order to fix the parameter b^. For 
example, the value b^ + 1 = 3.078 x 10"^ gives (Am^, sin^ 26) ~ (0.605 x 10"^ eV^ 0.9997), 
which is safely in the allowed region obtained by Faid et al. |^2| by using the solar neutrino 
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data from Super-Kamiokande and Borexino. In Table I, the neutrino oscillation parameters 
Amfj- and sin^ 26ff are listed. Here, the neutrino oscillation P(t'Q, — > Ujs) is given by 

P{u^ ^ .,(.|)) ^ E sin^ 2ef sin^ LAm^^H ^ ^ ^^^2 ^^a, ^-^2 , (45) 

where Am^(z/f'^) = |m^(z/f^) — m^(z/ff )|, Am^^ = m^(z/f'*) — rn^ivY), 



sin' 2ef = -4Re , (46) 



sin^ 29ff = -4Re 



and we have used the approximation m^(z/f_^) — m^(z/J^) ~ m^(z/f_^) — m^(z/Jl) ~ m^(z/fi) — 
"^^(^'J+) — m?{i'll) — 7Ti^(z/Jl) = Am^j- (i,j = e, /i,r). Also, since Am^^ ^ ^^fie^ 
Aml^ ~ Am^g, in Table I, we have denoted the sum sin^ 26*°^ + sin^ 26^"^ instead of sin^ 26^°^ 
and sin^ 2^"f as the effective mixing parameter sin^ 26 for Aml^. 

As seen in Table I, at Am' ~ 0.605 x 10"^° eV^ the oscillation i^e i^x is caused not only 
by i^e K (si^^ = 0.941) but also by — > i'^ (sin^ 26*1^ = 0.055). The small parameter 
+ 1 is insensitive to the numerical results except for Am^(z/f^). For example, by taking 
+ 1 = 2.8 X 10"^, we obtain Aiv?{vl'^) = 1.0 x 10~^ eV^, while the other numerical results. 
Am' and sin^ 29, are almost unchanged. The value sin^ 29^^ = 0.055 at Am' 10~^ eV' is 
too large compared with the small angle solution sin 26* ~ 7 X 10^3 in the MSW solutions 
T9| , |20[| , we cannot accommodate the present model to the MSW solutions. 



Thus, the model can explain the solar and atmospheric neutrino data by the large mixings 
Ue ^ vl and i/^, respectively. On the other hand, in the appearance experiments 

—>■ ^13 (tt, P = e, fi,T), the effective mixing parameters sin^ 29°f defined by (47) are highly 
suppressed because of the cancellation between z/f_^ and ufl . Therefore, in spite of a suitable 
value of Am', Am'^ ^ 1.4 eV^, the present model cannot explain the neutrino oscillation 

— > Ve experiment by the liquid scintillator neutrino detector (LSND) |23] at Los Alamos 



because our prediction is sin 26*1^ = 1.7 X 10"^ although Geiser has explained the LSND 
data by the oscillation v^. <-> in his model with three sterile neutrinos The effective 
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mixing angles sin^ 26"^^ = 3 x 10"^ at Am^ ~ 22 eV^ and sin^ 26"^^ = 0.0008 at Am^ ~ 1.4 
eV^ are too small to be detected by CHORUS and NOMAD |25| experiments at CERN. 



The effective electron neutrino mass (m^) for the case in Table I is also suppressed: 



i=l 



2.8 X 10"^ eV , (48) 



which is safely small compared with the upper bound on {niy) from the neutrino double j3 
decays {m^) < 0.68 eV |2| . 



The numerical results in Table I should not be taken rigidly. These values are sensitive 
to the parameter values used here, especially the value Am^(i/|"') is to the parameter (fej, + 
1). The values are only an example. We can give any values as far as these values are 
within the similar orders. Furthermore, we can also obtain similar results by using another 
parametrization 

m'^ = niLe"'^ , m'^ = niRe""^ , (49) 

instead of (30) (the numerical values of and ept are different from those in the parametriza- 
tion (30)). 

Generally, a model with sterile neutrinos is stringently constrained by the big bang 
nucleosynthesis (BBN) scenario. However, it has recently been pointed out by Foot and 
Volkas p3 that the stringent bounds can considerably be reduced by the lepton number 



asymmetry in the early Universe. Our pseudo-Dirac neutrinos z/^^ and i>^% have masses 
greater than 1 eV, and the squared mass difference Am^^ is Am^^ ~ 22 eV^. The mixing 
angles except those for Ui ^ z/f are sufficiently small. Therefore, the present model does not 
spoil the BBN scenario. 

In conclusion, we have discussed a possible neutrino-mass-generation scenario within the 
framework of the universal seesaw mechanism, especially, on the basis of a model with 
special matrix forms of rriL, Tn^ and Mp-, which can answer the question why only top 
quark t acquires the mass of the order of the electroweak scale = 0{mL), and can give 
reasonable quark masses and mixings in terms of the charged lepton masses. The model 
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provides three sets of the hght pseudo-Dirac neutrinos (z/f_^,t'ff) (i = e,fi,T) are mixing 
states uf^ ~ (z/j ± vl)l\pl between the conventional left-handed neutrinos Vi and sterile 
neutrinos z/f = [Nu ~ ^ri) I v^- The mixings z/g ^ and z/^ z/^^^ can give reasonable 
explanations of the solar and atmospheric neutrino data, respectively, with m(z/^^) ~ 5 
eV. Note that in spite of such a large mixing in a disappearance experiment, the effective 
mixing parameters sin 2^,f in the appearance experiments I'p are highly suppressed. 

This is the most remarkable feature in the model with the pseudo-Dirac neutrinos. If future 
appearance experiments z/^ at MINOS ||28|, ICARUS and K2K rule out the 



region (Am2,sin2 2^) ~ (lO^^ - IQ-^ eV^O.8 - 1) which is suggested by the atmospheric 
data, the present scenario with the three sterile neutrinos z/f will become an promising 
candidate of the possible interpretation of the discrepancy. 

At present, except for the LSND data, there is no positive motivation for such a model 



with three sterile neutrinos. Since the recent KARMEN data |^ seem to disfavor the 
LSND solution (Am^, sin^ 20) = (0.3 eV^ 0.04) - (2 eV^ 0.002), the motivation of the model 
with the sterile neutrinos seems to disappear. However, the present model with three sterile 
neutrinos is motivated not on the phenomenological interests in the neutrino data, but on 
the phenomenological success [§,0 of the universal seesaw mass matrix model for quarks. 
As far as the model is based on the SU(2)LxSU(2)/i;xU(l)y symmetry, the model includes 
the sterile fermions A^j, so that we inevitably have three sterile neutrinos vl in the model. If 
we do not want such the light sterile neutrinos z/f , we must consider an additional mechanism 
which guarantees m'^ = m'^ = 0. 

The present scenario can bring fruitful phenomenology into the neutrino physics, but 
there are many adjustable parameters in the present model, and those parameters have been 
determined by hand and by way of trial. In the present scenario, the active neutrinos z/j and 
sterile neutrinos vl are massless in the limit of A"^, *^ N'^ symmetry, and the small neutrino 
masses are induced by the violation of the A^^ ^ Nf^ symmetry, i.e., by the parameters El ^ 
and Er ^ 0. In other words, the neutrino masses are originated in "non-seesaw" mechanism 
in spite of the scenario in the "seesaw" model. We consider that the deviation between 
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m'^ {m'pi) and (^j?) can come from an asymmetric evolution of the Yukawa coupling 
constants from a unification scale fi = Ax to the electroweak scale /i = A^. However, 
whether this conjecture is reasonable or not is open question at present. 

Note added. After completion of this manuscript the authors became aware of a paper 
by Boweres and Volkas ( [hep-ph / 98043 1 0] ) , in which a neutrino mass matrix model with 
three pseudo-Dirac neutrinos has also been proposed within the framework of the universal 
seesaw mass matrix model. They have derived the light pseudo-Dirac neutrinos in their 
model with m'^ = m'^ = 0. Thier mass matrix resembles the model (5), but they assumed 
that 0{Mm) ^ O^Md). Even if we accept the assumption 0{Mm) ^ O^Md), their idea 
is not applicable to our model, because their light pseudo-Dirac neutrinos have masses of 
the order of m^M^^m/j in our scheme, so that the pseudo-Dirac neutrinos will acquire the 
masses of the same order as the charged leptons and light quarks have. However, it is 
noticeable that the model with m'^^ = m'p, = can also lead to three pseudo-Dirac neutrinos. 
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TABLES 

TABLE I. Neutrino oscillation parameters for the case: 6,^ + 1 = 3.078 x 10"^, 
= = 3.79 x IQ-^^, and A = 4.25 x 10^ {n = 8.50 x 10^) . 



Am2(i/f ) = 
0.605 X 10-1° eV^ 



4.97 X 10-=^ eV^ 



Am2(i/f ) = 
1.98 X 10-2 eV^ 



1.43 eV^ 



21.9 eV^ 



sin^ 261^ = 



sin' 29 f^^ 



sin^ 26^^ = 



sin^ 2^?f + 



sin^ 29"^ - 



0.941 

1.4 X 10-5 
2.7 X 10-9 



1.4 X 10-5 

0.940 
8.2 X 10-^ 



3.2 X 10-9 
6.7 X 10-^ 
0.992 



6.1 X 10-^ 

-3.4 X 10-9 
3.1 X 10-9 



-8.6 X 10-1° 

5.9 X 10-^ 
8.9 X 10-^ 



0.055 
-1.8 X 10-5 
0.0032 



0.00023 

0.0034 
-9.0 X 10-^ 



-2.1 X 10-1° 

0.00022 
9.1 X 10-^ 



1.7 X 10-^ 

0.00076 
3.6 X 10-6 



1.1 X 10-9 

3.0 X 10-6 
1.4 X 10-s 



0.9997 

0.0031 
1.0 X 10-5 



6.2 X 10-^ 

0.892 

1.3 X 10-5 



9.0 X 10-1^ 

5.1 X 10-s 
0.986 



0.00099 

0.210 
5.3 X 10-5 



3.8 X 10-6 
0.00091 
0.029 
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FIGURES 

FIG. 1. Mass spectrum of the pseudo-Dirac neutrinos uf^. The parameter values are the same 
as those in Table I. Boxes correspond to different mass eigenstates. The sizes of different regions 
in the boxes determine flavors of mass eigenstates, li/ajp. White and black regions correspond to 
the active and sterile flavors, respectively. 



ps 

Uq = 14+ 



ps 



4.84 eV 



Ami 



0.020 eV^ 



21.9 eV^ 



z/4 
1^1 



ps 



= V. 



ps 
e+ 




H 



1.20 eV T 0.005 eV^ 

1.43 eV^ 
0.078 eV T 6 X lO'^^ eV^ 
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